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It is shown that early suggested derivation of the Boltzmann kinetic equation for dilute hard 
sphere gas from the time-reversible BBGKY equations is incorrect since in fact a priori substitutes 
for them definite irreversible equations. Alternative approach to analysis of the hard sphere gas 
is formulated which conserves the reversibility and makes it clear that at any gas density one can 
reduce the BBGKY equations to the Boltzmann equation only in case of spatially uniform gas. 



PACS numbers: 05.20.-y, 05.40.-a 

1. The Boltzmann kinetic equation (BE) f]}-^ is 
one of most beautiful and fruitful models of theoretical 
physics. However, its status from viewpoint of statistical 
mechanics still stays under question. By the conventional 
opinion, solutions of BE coincide with solutions of exact 
equations of statistical mechanics, i.e. Bogolyubov-Born- 
Green-Kirkwood-Yvon (BBGKY) equations (BBGKYE) 

0, at least in the low-density gas limit or in the mathe- 
matically equivalent Boltzmann- Grad limit (BGL) when 
a— 7-0,z^— >-oo, /U = a^v — >■ , A = (tto^i/)"^ = const , 
with 1/ being mean gas density (concentration of gas par- 
ticles) and a and A being interaction radius and mean 
free path of gas particles, respectively [1, Q. The at- 
tempt to prove this assumption by considering the hard 
(elastic) sphere gas was made by Lanford and is known 
as "Lanford theorem" 6- 8] . But the Lanford result does 
not seem to be quite convincing because it was based on 
the formal series of iterations of BBGKYE which con- 
verges for absurdly small evolution times only, t < t 
(t ~ A/ \JT jm is mean free path time). Nevertheless, 
the "Lanford theorem" hardly is compatible with results 
of [gj (see also [IMl) where for the gas of "soft elas- 
tic spheres" it was shown that in case of its spatial non- 
uniformity BE does not follow from BBGKYE even under 
BGL. 

The aim of the present paper is to reveal the ori- 
gin of this contradiction starting from the "hard sphere 
BBGKY hierarchy" @, i, i, [H, [l3 . We will see that such 
the method for building solutions to this hierarchy as ap- 
plied in [B-Il] destroys its symmetry in respect to time 
inversion and insensibly replaces it by definite kinetic, 

1. e. irreversible, equations. Therefore the results of 

can not be qualified as BE derivation from BBGKYE. 
Besides, we will suggest and discuss a new approach to 
analysis of true solutions of the hard sphere BBGKY hi- 
erarchy. 

2. Let us consider the hard sphere gas @, [l^- There 
are no rigorous rules for transition to it from gas with 
smooth inter-particle interaction (thought a non-rigorous 
procedure was considered e.g. in [III). But it is possi- 



ble first to postulate the Liouville equation for the hard 
sphere system, as combination of Liouville equation for 
free particles, 

9F/at = -X;.v,-V,F at \v,-Vj\>a 

(where — Tpi/m , V,; — d/dvi ), and boundary condi- 
tions to it. 



F(...p*...p*...) =F(...p,...p,...) at r,- 

Pi,j ± VL (f7-(pj - Pi)) 



ail 



(1) 

( f2 is unit vector), which establishe continuity of (den- 
sity of) probability measure F along phase trajectories of 
particles under their collisions. Then, second, from here 
one can in usual way [H| deduce the desirable BBGKYE: 



dt 



^Vj- Vj Fn + i^^Ij^„+l Fn+1 , (2) 



where satisfaction of the boundary conditions ([T|) is pre- 
sumed, and the "collision operators" are defined by 



f^-(vfc-Vj) F{rk = Tj+aQ) dpk dn (3) 



It should be emphasized that these equations, 
like BBGKY in general, are time-reversible: if 
{Fn{t, r, p)} is some solution to equations (P)-® then 
{F„(— t, r, — p)} also is their solution. 

At this point the serious question does arise: how we 
have to deal with the conditions ([IJ? For the first look, 
we can merely use these conditions to express probabili- 
ties of post-collision ( out-) states via probabilities of pre- 
coUision ( in-) states and after that exclude the conditions 
from consideration. 

Then the collision operators take the form 

%k ^ = / Ji-^-i^k - V,)) 9{-n{vk - V,)) X 

x[F{rk = r,-an,p*,pl) ~ (4) 
- F{rk = Vj + afl , pj , pfe ) ] dpk dVL , 
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where 6Q is the Heavyside function indicating that the 
integration includes m-states only. 
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Just such transformed equations conventionally are as- 
sumed as a basis of the theory. All the more, they a priori 
are well predisposed to the Boltzmann's Stofizahlansatz 
[ij. Indeed, if we truncate the transformed hierarchy of 
equations at n — s , neglecting (s + l)-particle correla- 
tions, then the residuary s equations, - e.g. 

^ = -vi-Vi Fi + va^jdv2 Jdn (r!-V2i) 0(r!-V2i) x 

f)F ^ ^ f f 

-gf = -Yl ^r^3 ^2 + ^ J^a' / dps I dn (n-vsj) x 

X ein-V3j) [F2(p*)Fi(r3- r2i - an,pl) - (5) 
- F2{pj) Fi{r3= r2i +af7,p3)] 

at s = 2 [i^l , - under the BGL directly lead to the 

^ = -vrVi Fi + va^ j J (Q-w^i) e{n-V2i) x (6) 
X [ i^i (p^ ) i^i (p; ) - Fi (pi ) Fi (p2) ] dn dp2 



(we introduced 



— Vj ). This statement is in 



essence equal to the 'Lanford theorem" [29|_, because be- 
havior of any term of the iteration series [6|,|8| under BGL 
is determined by a finite number of the transformed equa- 
tions. 

3. Notice, however, that the transformation of EqslH 
by substituting (|4]) for ([3]), takes away from them their 
aforesaid time-reversibility property: now they are not 
invariant in respect to the inversion t p —> — p . 

Consequently solutions of the transformed equations in 
general are not solutions of the initial BBGKYE, and vise 
versa. In other words, that non-equivalent hierarchies of 
equations! 

Or, to be more precise, probably, we would be able 
to consider them as equivalent if we were able to deal si- 
multaneously with all possible many-particle correlations 
(of arbitrary high orders). But they are certainly non- 
equivalent if we neglect even a part of correlations [s^ . 

Therefore the term by term consideration of iteration 
series for the transformed equations [6, 8] leads away from 
a true solution of the initial "hard sphere BBGKY hier- 
archy" ([I])-©. Instead it offers solution of another, in 
essence, postulated kinetic equations. Seemingly, this re- 
sult sorts with true BBGKY solutions approximately like 
complex-energy solutions of the Schrodinger equations 
(or, in classical statmechanics, solutions of the Li- 
ouville equations corresponding to the complex "Ruelle- 
PoUicott resonances" 1191) sort with their actual solu- 



tions [31 1 . 

4. Thus, the replacement of ([3]) by gives the same 
effect as the Boltzmann's Stofizahlansatz. From the phys- 
ical point of view, this "ansatz" is bad as for it forces us 
to neglect pre-collision (m-) inter-particle correlations 
and hence neglect fluctuations in "relative frequency" , 
or "frequency ratio" ("chastost"' in Russian f^Tl) of col- 
lisions 



Indeed, any collision starts from a pre-collision config- 
uration and finishes with a post-collision configuration, 
hence, fluctuations in "relative frequency of collisions" 
equally give rise to both out- and m -correlations [s?] . 

It is necessary to emphasize that we say about statis- 
tical correlations which do not presume presence of some 
cause-and-consequence relations beyond them fssj . 

The cause of the "fluctuations in relative frequency 
of collisions" is mere absence of back reaction to them 
when they do not disturb the system's state (for instance, 
when relative frequencies of mutually time-reversed colli- 
sions fluctuate with keeping deflnite proportions between 
them) [sl] . 

Clearly, these fluctuations are as well reflected by the 
distribution functions (DF) {Fn(t}} as strong is spatial 
non- uniformity of the system d, [l^] , and therefore they 
are reflected in the form of correlations between particles' 
coordinates (while their velocities can be uncorrelated as 
in the Boltzmann's theory) [35| . 

These spatial correlations [36| , in their turn, do mean 
that the DF values at collision conflgurations, e.g. 
F2(r2 = ri + ail , represent independent on Fi{t) and 
complementary to Fi (t) characteristics of statistical en- 
semble ^ . 

5. In view of the aforesaid, we have to come back to 
the question how we must deal with the conditions 

Since "contact" DF's values which enter ([1]), first of all 
F2(ri — r2+ ail , play the role of "governing parameters" 
for BBGKY hierarchy ^ as the whole, it is natural to 
treat them as independent on Fi characteristics of gas. 
In more detail, when considering F2(ri — r2 + ail, we 
inevitably come to rest (in the collision integral) against 
three-particle configurations corresponding to pairs of in- 
finitely close [air collisions. The, considering such config- 
urations, we will come to analogous four-particle ones, 
and so on. Categorizing all them, one would construct a 
full (infinite) system of equations for the "contact" Dfs. 
Such a system, of course, would be time-reversible fssi] . 

Realization of such a program just would give the an- 
swer to the question. On this way, one can easy im- 
mediately see a mechanism of generation of the spatial 
correlations and destroying the Stofizahlansatz. 

Making the first step, let us rewrite the second of 
BBGKYE in the "pseudo-Liouville" form Q: 

BF f 

= a^ (vi2 -l^) (5(ri2 - ail) F^ dil ~ 



-Vl2- 



dF2 
dri2 



V1+V2 dF2 



(7) 



Here, the coordinates r^- may enter the forbidden region 
|ri — r2| < a, where F2 = , and the new (first on 
r.h.s.) term represents a force of repulsion of particles at 
the border of this region [s^ . Besides, we separated the 
relative displacement of particles, with ri2 = ri — r2 , 
and motion of their center of mass, R = (ri -f r2)/2. 
The first of these two in turn can be divided into norma; 
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and tangential components: 



dF2 , „, dF^ 



dr 



12 



Opi2 \Pl2 



where = ri2/|ri2| , pi2 = |ri2| and 



A(r2)f = f - r2(Sl-f) 



[f7 X [r2 X f ]] 



Next, consider, with the help of the contact DF 

F2(ri = r2 + ail) = F2''\t, R, fi, pi, P2) , taking in mind 
that the first term on r.h.s. of ([7]) and the first (nor- 
mal) component of ([5]) have quite similar, but oppositely 
signed, singularities which compensate one another. Let 
us assume that this is exact compensation, that is 



a^J (vi2f7') S{pi2n - aft') F2 dfl' - (vi2f7) 



dpi2 



(9) 



Then from ©-(ini) the necessary autonomous evolution 
equation for the pair contact DF does follow: 



dF 



V1+V2 dF. 



dK 



(10) 



J , 3 ^"3 



(c) 

where Fg is the mentioned contact DF for two "bound 
together" pair collisions. 

Notice that formally © is not assumption but iden- 
tity. It expresses continuity f probability distribution at 
breaks of phase trajectories because of collisions, that is 
the same as the condition ([T]) does express. In essence, 
this is analogue of equalities (3)- (4) from ^ for a gas 
with smooth interaction 1401 . 



From the equation ([TO]) it it clear that any spa- 
tial inhomogeneity, which induces dF^^^ /dvi ^ and 
(8) (9^2 '^''/9R ^ , automatically excludes possibility of re- 
duction of ^2^^^^ to Fi and thus BBGKYE to BE, abso- 
lutely independently on value of p, [4l| . 

Notice also that, firstly, the equation (|10l) is reversible 
and besides invariant in respect to replacing pi,P2 by 
p*,P2 , as it should be according to ^ (thus, the func- 
tion of ([!]) now is extension of this symmetry property 
from equations to their solutions). 

Secondly, the equality ([9]) says, in particular, that 
dF2/dpi2 =0 at pi2 = a-hO and (viafJ) ^ 0. This 
is natural analogue of behavior of gas density nearby a 
flipping (in accordance with ([1])) surface. 



6. The contact DF F2 serves as a measure of mean 
(ensemble averaged) number density of pair collisions. It 
is clear that it drifts with the center of mass velocity of 

ic] 

colliding particles. Similarly, DFs Fs will drift with 
velocities (vi -I- ... -I- Vs)/s. Therefore they are mutu- 
ally independent. All together they determine statistics 
of key s-particle configurations which produce all other 
configurations and eventually evolution of -Fi [43 . 

7. Formulation of equations for Fg'^' , F4 , etc. 
we leave for the future. At present, it is more impor- 
tant to point out existence of alternative, just discussed, 
treatment of the hard sphere BBGKY hierarchy. This 
treatment, in contrast with the conventional one, does 
not ignore the fundamental reversibility property of the 
BBGKY equations but uses it as basis of definite con- 
structive approach to solutions of these equations. A 
qualitatively similar approach to a gas of "soft spheres" 
was tested in citeil,i2,pl2. The model of "hard spheres" 
is of special interest because, expectedly, it can more eas- 
ily achieve formal quantitative rigor. 
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in application to other systems in in [ly, [2J, [2g], and 
besides recently in [l^ - [T^ . In it was shown that these 
fluctuations are indifferent to a degree of smallness of 
the gas parameter fi . In opposite, just at ^ (in 
BGL) the "absence of back reaction" is especially easy 
understandable. 

There is a key to understanding the 1/f-noise observed 
in various physical systems [gHlll. |23 - [25| ]. And we can 
expect that exact time-reversible solutions of BBGKYE 
contain 1/f fluctuations in kinetic characteristics of the 
system. 
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[37] Already because for ensemble averages generally the in- 
equalities {u") 7^ {v)" take place, where p is local gas 
density. 

[38] Just by this reason it would present true statistical 
weights of any of kinematically possible scenarios of col- 
lisions. 

Generally, since irreversibility equally manifests itself in 
both opposite time directions, its completely adequate 
description can be done only by reversible equations! 

[39] In place of (Vi2$(ri2)) • (9^2 /Spi - 9^2 /(9p2) in case 
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surprising: regardless of velocity of a given particle, in its 
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[42] In more detail, on such kind of DFs see [9l-[Tll]. To avoid 
misunderstandings, it is useful to underline that the 
smallness of probabilities of contacts of two or more par- 
ticles in no way means smallness of corresponding contact 
DFs since the latter (as well as all the DFs under consid- 
eration) represent density of probability. 



